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Abstract 

The production of charm quarks is studied in deep-inelastic electron-photon scattering using 
data recorded by the OPAL detector at LEP at nominal e + e~ centre-of-mass energies from 
183 to 209 GeV. The charm quarks have been identified by full reconstruction of charged D* 
mesons using their decays into D°7r with the D° observed in two decay modes with charged 
particle final states, Kir and Ktttttt. The cross-section cr D for production of charged D* in 
the reaction e + e~ — > e + e~ D* X is measured in a restricted kinematical region using two 
bins in Bjorken x, 0.0014 < x < 0.1 and 0.1 < x < 0.87. From cr D the charm production 
cross-section cr(e + e~ — > e + e~ ccX) and the charm structure function of the photon F^ c are 
determined in the region 0.0014 < x < 0.87 and 5 < Q 2 < 100 GeV 2 . For x > 0.1 the 
perturbative QCD calculation at next-to-leading order agrees perfectly with the measured 
cross-section. For x < 0.1 the measured cross-section is 43.8 ± 14.3 ± 6.3 ± 2.8 pb with a 
next-to-leading order prediction of 17.0 1|J pb. 
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1 Introduction 



Much of the present knowledge of the structure of the photon has been obtained from deep- 
inelastic electron-photonQ scattering at e + e~ colliders 0. With the high statistics available 
at LEP2, it is possible to investigate the flavour composition of the hadronic structure 
function K7. The easiest flavour component of F% to measure directly is the charm part 
F2 C , because charmed hadrons like the D* meson are produced with a large cross-section 
and can be identified by the well established method of full reconstruction of decays of 
charged D* mesons. 

The measurement is based on deep-inelastic electron-photon scattering, e(k) 7(p) — > 
e(k') ccX, proceeding via the exchange of a virtual photon, 7* (9), (the symbols in brackets 
denote the four- momentum vectors of the particles). Experimentally, for these events one 
electron is observed in the detector together with a hadronic final state, and the second 
electron, which is only slightly deflected, leaves the detector unobserved. The determina- 
tion of F^c exploits the fact that the differential cross-section of this reaction as function of 
Q 2 = —q 2 and Bjorken x, defined as x = Q 2 /(2p ■ q), is proportional to F^x, Q 2 ). 

For a given gluon distribution function, the contribution to F% from charm quarks can 
be calculated in perturbative QCD, thanks to the sufficiently large scale established by the 
charm mass, and predictions have been evaluated at next-to-leading order (NLO) accuracy 
in 0. As for light quarks, K7 C receives contributions from the point-like and the hadron-like 
components of the photon structure, as shown schematically in Figure |I[ These two contri- 

1 If not mentioned explicitly charge conjugation is implied, and for conciseness positrons are also referred 
to as electrons. Natural units h — c — I are used throughout. 
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Figure 1: Examples of leading order diagrams contributing to a) the point-like, and b) the 
hadron-like part of F^. 



butions are predicted to have different dependences on x, with the hadron-like component 
dominating at very low values of x and the point-like part accounting for most of at 
x > 0.1. 

The analysis presented here is an extension of the measurement of F^ presented in ||, 
using basically the same analysis strategy, but a much larger data sample and refined Monte 
Carlo models. It is based on data recorded by the OPAL experiment in the years 1997-2000, 
with an integrated luminosity of C = 654.1 pb _1 for nominal e + e~ centre-of-mass energies, 
y/s^, from 183 to 209 GeV, with a luminosity weighted average of \fs^ = 196.6 GeV. 



2 The OPAL detector 



A detailed description of the OPAL detector can be found in |4j], and therefore only a brief 
account of the main features relevant to the present analysis is given here. 

The central tracking system is located inside a solenoidal magnet which provides a uni- 
form axial magnetic field of 0.435 T along the beam axisQ The magnet is surrounded 
by a lead-glass electromagnetic calorimeter (ECAL) and a hadronic sampling calorimeter 
(HCAL). Outside the HCAL, the detector is surrounded by muon chambers. There are sim- 
ilar layers of detectors in the endcaps. The region around the beam pipe on both sides of 
the detector is covered by the forward calorimeters and the silicon-tungsten luminometers. 

Starting with the innermost components, the tracking system consists of a high precision 
silicon microvertex detector ||, a precision vertex drift chamber, a large volume jet chamber 
with 159 layers of axial anode wires, and a set of z chambers measuring the track coordinates 
along the beam direction. The transverse momenta pt of tracks with respect to the z direction 

2 In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis points 
upwards and the z axis points in the direction of the electron beam. The polar angle 6 and the pseudorapidity 
T) = — m(tan(0/2)) are defined with respect to the z axis. 
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of the detector are measured with a precision of cr Pt /p t = y 0.02 2 + (0.0015 • p t ) 2 (pt in GeV) 
in the central region, where |cos#| < 0.73. The jet chamber also provides energy loss, 
dE/dx, measurements which are used for particle identification. 

The ECAL covers the complete azimuthal range for polar angles satisfying | cos#| < 0.98. 
The barrel section, which covers the polar angle range | cos^l < 0.82, consists of a cylindrical 
array of 9440 lead-glass blocks with a depth of 24.6 radiation lengths. The endcap sections 
(EE) consist of 1132 lead-glass blocks at each end with a depth of more than 22 radiation 
lengths, covering the polar angle region of 0.82 < | cos#| < 0.98. 

The forward calorimeters (FD) at each end of the OPAL detector consist of cylindrical 
lead-scintillator calorimeters with a depth of 24 radiation lengths divided azimuthally into 16 
segments. The electromagnetic energy resolution is about 18%/y/E (E in GeV). The accep- 
tance of the forward calorimeters covers the angular range between 47 and 140 mrad from the 
beam direction. Three planes of proportional tube chambers at 4 radiation lengths depth in 
the calorimeter measure the direction of electron showers with a precision of approximately 
1 mrad. 

The silicon tungsten detectors (SW) |j) are located in front of the forward calorimeters at 
each end of the OPAL detector. Their clear acceptance covers a polar angular region between 
33 and 59 mrad. Each calorimeter consists of 19 layers of silicon detectors and 18 layers of 
tungsten, corresponding to a total of 22 radiation lengths. Each silicon layer consists of 16 
wedge-shaped silicon detectors. The electromagnetic energy resolution is about 25%/\/E (E 
in GeV). The radial position of electron showers in the SW calorimeter can be determined 
with a typical resolution of 0.06 mrad in the polar angle 6. The SW detector provides the 
luminosity measurement. 



3 Monte Carlo simulation 

The Monte Carlo models HERWIG6.1 and PYTHIA6.1 |§, both based on leading order 
(LO) matrix elements and parton showers, are used to model the deep-inelastic electron- 
photon scattering events e + e~ — ► e + e~cc X. For both Monte Carlo models, the charm quark 
mass is chosen to be m c = 1.5 GeV. 

In HERWIG6.1, the cross-section is evaluated separately for the point-like and hadron-like 
contributions using matrix elements for massive quarks, together with the GRV parametri- 
sation |J for the gluon distribution of the photonf]. This is an improvement compared to 
the analysis presented in ||, where HERWIG5.9 was used with matrix elements for massless 
quarks only. For that model the effect of the charm quark mass was accounted for only rather 
crudely by not simulating events with hadronic masses of less than 2m c . The fragmentation 

3 The massive matrix elements for the point-like contribution have been implemented into the standard 
HERWIG6.1 model by J. Chyla. The hadron-like component is based on the matrix elements of the boson- 
gluon fusion process. 
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into hadrons is based on the cluster model for HERWIG6.1, using the OPAL tuneQ. 



To obtain an estimate of the dependence of the measurements on the details of the mod- 
elling of D* production and fragmentation, a second model, based on the PYTHIA program, 



is used. The spectrum of photons with varying virtualities is generated in PYTHIA [10 
The point-like and the hadron-like processes, which in PYTHIA are denoted by direct and 
resolved, are then simulated separately using the matrix elements for the production of mas- 
sive charm quarks, 77 — > cc (subprocess ISUB=85) and 7g — > cc (subprocess ISUB=84). 
Since these matrix elements are valid only for real photons (Q 2 = 0) the Q 2 dependence 
of the cross-sections is not expected to be correctly modelled. The cross-sections for both 
contributions are therefore taken from HERWIG, and PYTHIA is only used for shape com- 
parisons. In PYTHIA the fragmentation into hadrons is simulated using the Lund string 
model. 

All Monte Carlo events were passed through the full simulation of the OPAL detector [ |TT|j . 
They are analysed using the same reconstruction algorithms as applied to the data. 



4 Kinematics and data selection 



To measure -F^ c (x, Q 2 ), the distribution of events in x and Q 2 is needed. These variables are 
related to the experimentally measurable quantities E tSug and # tag by: 

Q 2 = 2£ b £ tag (l-cosfl tag ) (1) 

and 



Q 2 + W 2 + P 2 ' w 

where is the energy of the beam electrons, E tag and 9 ta g are the energy and polar angle of 
the deeply inelastically scattered electron, W 2 is the invariant mass squared of the hadronic 
final state, and P 2 = —p 2 is the negative value of the virtuality squared of the quasi-real 
photon. It is required that the associated electron is not seen in the detector. This ensures 
that P 2 < Q 2 + W 2 , therefore P 2 is neglected when calculating x from Equation |2[ The 
electron mass is neglected throughout. 

Deep-inelastic electron-photon scattering candidate events are selected as follows: 

1. The calorimeter cluster with the highest energy in either SW or FD is taken as the 
electron candidate. The polar angle 8 ta , g is measured with respect to the original beam 
direction. An electron candidate is required to have E iag > 0.5E h and 33 < 6 tag < 
55 mrad (SW) or 60 < tag < 120 mrad (FD). 

4 The main changes are that meson states that do not belong to the L = 0, 1 supermultiplets are removed, 
and that the parameters CLSMR(l), PSPLT(2) and DECWT have been changed from their default values of 
0.0, 1.0, and 1.0 to 0.4, 0.33, and 0.7. A detailed description of the tune can be obtained from the HERWIG 
web interface. 



6 



2. To ensure that the virtuality of the quasi- real photon is small, the highest energy 
electromagnetic cluster in the SW and FD detectors in the hemisphere opposite the 
scattered electron must have an energy E a < 0.25 (the anti-tag condition). 

3. At least three tracks must be found in the tracking chambers. A track is required to 
have a minimum transverse momentum with respect to the beam axis of 0.12 GeV, 
and to fulfill standard quality criteria ||. 

4. To reduce background from e + e~ annihilation events with charged D* mesons in the 
final state, the sum of the energy of all calorimeter clusters in the ECAL is required 
to be less than 0.5E^,. Electromagnetic calorimeter clusters have to pass an energy 
threshold of 0.1 GeV for the barrel section and 0.25 GeV for the endcap sections. 

5. To reduce the e + e~ annihilation background further, the visible invariant mass of the 
event, W v i S , should be less than 0.65-Eb- The invariant mass W v - ls is calculated using 
the momenta of tracks and using the energies and positions of clusters measured in 
the ECAL, the HCAL, the FD and the SW calorimeters. Clusters in the SW and FD 
in the hemisphere of the tagged electron are excluded. A matching algorithm [12|] is 
applied to avoid double counting the energy of particles that produce both tracks and 
clusters. All tracks, except for the kaon candidate identified in the D* reconstruction, 
are assumed to be pions. 

6. A fully reconstructed charged D* candidate has to be present. The transverse momen- 
tum of the D* with respect to the beam axis, pf*, should fulfill pf* > 1(3) GeV for 
SW(FD)-tagged events, and the pseudorapidity r/ D * should be in the range |?7 D *| < 1.5. 
The D* mesons are identified using their decay into D°7r with the D° observed in two 
decay modes with charged particle final states, Ktc and Kimn. The following quality 
criteria are applied to the D° candidate. For the kaon candidate, the dE/dx proba- 
bility for the kaon hypothesis should exceed 10%, and for all pion candidates the pion 
probability should be above 0.5%. In addition, the mass of the D° candidate should 
lie in the window 1.79-1.94 GeV, and the cosine of the angle between the kaon mo- 
mentum in the D° rest frame and the D° momentum in the laboratory frame should 
be below 0.9. The cut on the cosine of the angle reduces the background which peaks 
at 1, whereas the pseudo-scalar D° decays isotropically in its rest frame. 

With these cuts 1653 events are selected in the data having Am < 0.21 GeV, where Am is 
the difference between the D* and D° candidate masses. They are shown in Figure |^. The 
number of events reconstructed in the signal region defined by 0.1424 < Am < 0.1484 GeV 
is 115. No event with more than one D* candidate in the signal region has been found in 
the data. 



Using Monte Carlo simulations, the expected background to the deep-inelastic electron- 
photon scattering events from all other Standard Model physics processes that potentially 
contain final state D* mesons is found to be about one event and is neglected. The back- 
ground from random coincidences between off-momentumQ beam electrons faking a scat- 
tered electron and photon-photon scattering events without an observed electron has been 
estimated to be below 1.4 events and has been neglected. Thus, only the combinatorial 

5 Off- momentum electrons originate from beam gas interactions far from the OPAL interaction region and 
are deflected into the detector by the focusing quadrupoles. 
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background from deep-inelastic electron-photon scattering events e + e~ — > e + e~qq with 
q = u, d, s and c has to be subtracted from the data. The production of bottom quarks 
is suppressed compared to charm production due to the larger mass and smaller charge and 
is negligible 



5 Results 



5.1 Comparison of data and theoretical predictions 

Figure shows the difference between the D* and D° candidate masses for both decay 
channels combined. A clear peak is observed around the mass difference between the D* and 
the D° meson, which is Am = 0.1454 GeV |I3|. The number of signal events, N® c , has 
been obtained from an unbinned maximum likelihood fit to this distribution. Similar fits are 
performed in two regions of x vis , calculated from Equation || using W vis and the measured 
value of Q 2 . The fit function contains a Gaussian for the signal and a power law function 
of the form a ■ (Am — m 7T ) b , for the background contribution. Because the number of signal 
events is small, the mean value of the Gaussian is kept fixed to the PDG value in the fit. The 
width and the normalisation are determined by the fit, which gives 55.3 ± 11.0 signal events, 
with 23.6±7.4 and 31.4±8.1 events for x v - ls < 0.1 and x v - ls > 0.1, where the uncertainties are 
statistical. The fitted width is 0.70 ± 0.16 MeV when using the whole sample. The quality 
of the fit is satisfactory. The x 2 between the fitted curve and the data in Figure |2| is 94 for 
71 non-empty bins. The corresponding y 2 for events with x v - ls < 0.1 and x v - ls > 0.1 are 71 
and 93. 

The result of the fit for all events is shown in Figure |2] together with the absolute pre- 
diction of the combinatorial background measured from the data using events with a wrong- 
charge pion for the D* — > D°7r decays. The fit agrees with this second estimate of the 
combinatorial background. The numbers of signal events for x v - ls < 0.1 and x v - ls > 0.1 pre- 
dicted by the HERWIG model are 12.1 ± 0.6 and 30.0 ± 1.1, where the uncertainties are 
statistical. 

Figure ^| shows, for the signal events, the distributions of four global event quantities, Q 2 , 
W v i s , x vis and the charged multiplicity iV ch , and two variables related to the kinematics of 
the D* candidates, p^* and |r/ D *|. The signal events are shown for the Am region 0.1424 < 
Am < 0.1484 GeV, after subtracting the combinatorial background in that region. The 
normalisation of the background is given by the result of the fit for the complete sample. 
The shape of the background distributions is taken from the data using both the events with 
the wrong-charge combinations and the events with the correct charge combinations that 
fulfill 0.16 < Am < 0.21 GeV, see Figure [2| It has been verified from Monte Carlo that 
in this range the distributions shown are independent of Am. Subtracting the background 
this way is statistically more precise than using only the events with the wrong-charge 
combinations in the signal region. The data are compared to the prediction of the HERWIG 
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and PYTHIA Monte Carlo models, normalised to the data luminosity. For the HERWIG 
prediction the hadron-like component (HERWIG HL) is shown on top of the point-like 
component (HERWIG PL). Overall, the HERWIG model describes the data distributions 
quite well, and the shape of the PYTHIA prediction is also consistent with the data. 

In what follows, the cross-sections has been calculated based on the HERWIG predic- 
tion, while the PYTHIA program was used as a second Monte Carlo model to estimate the 
uncertainty stemming from the Monte Carlo description of the data. 



5.2 Determination of the charm production cross-section 

The inclusive cross-section cr D has been extracted in a well-measured kinematic region where 
p^ > 1 GeV for an electron scattering angle of 33 — 55 mrad, or p® > 3 GeV for 60 — 
120 mrad, |?7 D *| < 1.5 and 5 < Q 2 < 100 GeV 2 , using almost the whole accessible Q 2 range 
defined by 9 ta , g and E t3jg . 

As in the previous investigation ||, the analysis is performed in two bins in x, 0.0014 < 
x < 0.1 and 0.1 < x < 0.87. The x range is limited by the Q 2 range, by the minimum 
kinematically allowed invariant mass W = 3.88 GeV needed for the production of a D* 
meson together with a second charmed hadron, and by the event selection cut on W v - ls . To 
take into account the detector acceptance and resolution in x the data are corrected using a 
2x2 matrix, which contains the information of the correlation between the measured value 
x vis and the generated x for the two bins in x as given by Monte Carlo. The effects due to 
the migration in Q 2 are small compared to those in x, and are neglected. 

For a given D* decay channel the selection efficiency for x < 0.1 is given by the ratio of 
the number of reconstructed D* mesons originating from events with x < 0.1 to all generated 
D* mesons in that channel with x < 0.1 in the restricted kinematic region defined above. 
For HERWIG6.1, the selection efficiencies for x < 0.1 are (16 ± 1)% and (12 ± 1)% for 
the point-like and hadron-like components. For x > 0.1 they are (19 ± 1)% and (19 ± 5)% 
respectively. Combining the two components yields efficiencies of (15 ± 1)% and (19 ± 1)% 
for x < 0.1 and x > 0.1. The corresponding selection efficiencies for the PYTHIA sample 
are (14 ± 1)% and (20 ± 1)%. This means that in each bin the predicted efficiencies from 
the HERWIG and PYTHIA models are consistent, so the size of their error of about 10% is 
taken as the systematic uncertainty on the efficiency. 

Table |T|a) summarises the measured values of a D for D* production in the restricted 
kinematic region defined above: 

_D* cor /o\ 

° = BR-£- (3) 

It is calculated from the number of D* events, A"°*, obtained from N^* c determined by 
the likelihood fit to the full sample, together with the 2x2 matrix correction based on 
the HERWIG model. The statistical correlation between the two data points is —28%. 
For the combined D* branching ratios for the 3-prong and 5-prong decay modes, BR = 
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7.66 ± 0.22 % Pi is used. 



The present measurement makes use of the estimation of the systematic error in || with 
some small changes. The statistical error from the background subtraction is treated differ- 
ently and is now included in the statistical error of the number of signal events determined 
by the fit. For each of the x regions, a total correlated systematic uncertainty of 10%, as 
evaluated in 0, is attributed to the sum of the uncertainties stemming from the branch- 
ing ratios, the imperfectness of the modelling of the central tracking detectors, the kaon 
identification and the variation of the energy scale of the electron candidate. 

For the determination of the cross-section of charm production, a(e + e~ e + e~ccX), 
the Monte Carlo models are used for extrapolation. This allows cr(e + e~ — > e + e~ cc X) to be 
calculated via the relation 

a(e + e" -> e + e" cc X) = — • a(e + e" -> e+e" D* X) 

v ; 2-/(c->D*) V ' 



1 



K MC -a B \ (4) 



2 • f(c - D* 



The value used for the charm to D* hadronisation fraction, f(c — > D*), is taken to be 
0.235 ± 0.007 ± 0.007, independent of x. It has been estimated in |TJ| by averaging results 
obtained in e + e~ annihilation at LEP1 as well as at lower e + e~ centre-of-mass energies. The 
extrapolation factor TZ MC is defined as the ratio of the number of all generated D* mesons 
in the full kinematic region of and |^ D | divided by the number of generated D* mesons 
in the restricted region. For x < 0.1 the TZ MC values for the HERWIG and PYTHIA models 
are 6.15 ± 0.03 and 6.23 ± 0.03, while for x > 0.1 the corresponding numbers are 4.73 ± 0.02 
and 4.59 ±0.02, where the errors are statistical. The central values of the extrapolated cross- 
sections are obtained using the HERWIG numbers, and the differences between the HERWIG 
and PYTHIA predictions have been taken as estimates of the extrapolation errors. 

In Figure f|a) and Table |l]b), the measured cross- sect ions <r(e + e~ — > e + e~ccX) are 
compared to the calculation of @ performed in LO and NLO and to the prediction from 
HERWIG. The NLO prediction uses m c = 1.5 GeV. The renormalisation and factorisation 
scales are chosen to be p% = = Q 2 . The calculation is obtained for the sum of the point- 
like and hadron-like contributions to F^, where for the calculation of the hadron-like part 
the GRV-NLO parametrisation is used. These GRV-NLO parton distributions of the photon 
have been found to describe the OPAL jet data [TE\. The NLO corrections are predicted 



to be small for the whole x range. The NLO calculation is shown in Figure [| as a band 
representing the uncertainty of the theoretical prediction, evaluated by varying the charm 
quark mass between 1.3 and 1.7 GeV and by changing the renormalisation and factorisation 
scales in the range Q 2 /4 < [i 2 = fi 2 < 4Q 2 , taking the largest difference from the central 
value as the error. The HERWIG prediction falls within the uncertainty band of the NLO 
calculation. 



For x > 0.1, the predicted NLO cross-section in Figure |]a) agrees well with the data. 
For x < 0.1, the situation is different. The NLO calculation predicts the hadron-like and 
point-like component to be of about equal size, and the sum is smaller than what is observed 
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for the data. 

The point-like contribution can be calculated with small uncertainties, i.e. in LO it de- 
pends only on m c . The prediction for the hadron-like part is more uncertain, especially 
because it depends on the gluon distribution of the photon, for which experimental infor- 
mation is limited. Given the good agreement at large x between the data and the NLO 
prediction, dominated by the point-like component, the NLO point-like prediction has been 
subtracted from the measured cross-section in the region 0.0014 < x < 0.1. This leads to 
the value of 34.5 ± 14.3 ± 6.9 pb for the hadron-like contribution to the cross-section, which 
is to be compared with the NLO prediction of 7.7 pb, when using the GRV-NLO parton 
distributions of the photon. 



5.3 Extraction of F^ c 

The value of the charm structure function F2 tC (x,Q 2 )/a of the photon, averaged over the 
corresponding bin in x, and given at a fixed value of Q 2 , is determined by 

FUx,Q 2 )/a = a(e+e- eVccI) ■ ( - f/ X > - ) , (5) 



cx(e + e~ — > e+e~ ccX) 



NLO 



where the ratio (F^ c (x, Q 2 )/a/a(e + e~ — ► e + e~ccX)) NLO is given by the NLO calculation 
of . This approach assumes that for a bin of x the ratio of the structure function, averaged 
over x and evaluated at the average Q 2 value of the data, and the cross-section within the 
region of Q 2 , are the same for the data and the NLO calculation. The measurement is 
given at Q 2 = 20 GeV 2 , which roughly corresponds to the average value of Q 2 observed for 
the data. The F2 :C (x,Q 2 )/a values are listed in Table |I|c) and shown in Figure f|b) on a 
logarithmic scale in x. 

In addition to the structure function of the full NLO calculation, the predicted hadron- 
like component is also shown in Figure |]b). This contribution is very small for x > 0.1 and 
therefore in this region the NLO calculation is a perturbative prediction which depends only 
on the charm quark mass and the strong coupling constant. This prediction agrees perfectly 
with the data. 



To illustrate the shape of F^ c the data are also compared to the GRS-LO [HI prediction 
and to the point-like component alone, both shown for Q 2 = 20 GeV 2 . The point-like 
contribution is small at low x. The data span a rather large range in Q 2 , in which the 
change of the predicted K7 C is large. The maximum value of _F£ C for x > 0.1 rises by about 
a factor of five between the lower boundary of Q 2 = 5 GeV 2 and the higher boundary of 
Q 2 = 100 GeV 2 . In addition, the mass threshold, W = 2m c , introduces a Q 2 dependent 
upper limit in x which, in the Q 2 range studied, varies from x = 0.35 to x = 0.9. 

To be able to compare the data directly with the curves from the GRS-LO predictions the 
points are placed at those x positions that correspond to the average predicted F^ c within the 
bin [O . These x values are calculated both for the full F^ c and for the point-like component 
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alone. The data points are located at the mean of the two values and the horizontal error 
bar indicates their difference. For x > 0.1 the difference is invisible in Figure |]b). 

For x < 0.1 the hadron-like contribution to F^ c (x,Q 2 )/a, determined as for the cross- 
section, amounts to 0.154 ± 0.059 ± 0.029, to be compared with the NLO prediction of 
nn 2 fi+ 0-007 

U-UZO _ Q.005- 



6 Conclusions 



The production of charm quarks is studied in deep- inelastic electron-photon scattering using 
data recorded by the OPAL detector at LEP at nominal e + e~ centre-of-mass energies from 
183 to 209 GeV in the years 1997-2000. The result is based on 654.1 pb _1 of data with 
a luminosity weighted average centre-of-mass energy of 196.6 GeV. The measurement is an 
extension of the result presented in H, using basically the same analysis strategy, but with 
improved Monte Carlo models and higher statistics. The two OPAL results are consistent 
and the new measurement supersedes the result in [|J. 

The cross-section <r D * for D* production in the reaction e + e~ — > e + e~ D* X is measured 
in the deep-inelastic scattering regime for a restricted region where p^* > 1 GeV for an 
electron scattering angle of 33 — 55 mrad, or > 3 GeV for 60 — 120 mrad, |?7 D *| < 1.5 
and 5 < Q 2 < 100 GeV 2 , divided into two bins in x, 0.0014 < x < 0.1 and 0.1 < x < 0.87. 
Within errors the cross-sections can be described by the HERWIG model. From cr D the 
cross-section cr(e + e~ — > e + e~ cc X) is obtained by extrapolation in the same bins of x and 
Q 2 , using Monte Carlo. 

The charm structure function F^^x, Q 2 )/a of the photon is evaluated at Q 2 = 20 GeV 2 
for the same bins of x. For x > 0.1, the perturbative QCD calculation at next-to- leading 
order agrees perfectly with the measurement. For x < 0.1 the point-like component, however 
calculated, lies below the data. Subtracting the NLO point-like prediction, a measured value 
for the hadron-like part of 0.154 ± 0.059 ± 0.029 is obtained. 
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a) 


a D * [pb] 


0.0014 < x < 0.1 


O.K x < 0.87 


OPAL 
HERWIG 


3.1 ± 1.0 ±0.5 
1.13 (0.71 + 0.43) 


2.6 ±0.9 ±0.3 
2.05 (2.02 + 0.03) 



b) 


cr(e + e — > e + e ccX) [pb] 


0.0014 < x < 0.1 


O.K x < 0.87 


OPAL 
HERWIG 
LO 
NLO 


43.8 ± 14.3 ±6.3 ±2.8 
17.0 (8.5 + 8.4) 
16.0 (8.1 + 7.9) 

i7.oil;l (9.318:1 + 7.71 ?;§) 


26.2 ± 8.8 ± 3.2 ± 1.3 
27.4 (27.1 + 0.4) 
26.7 (26.3 + 0.4) 
30.8 1^ (30.4±«:i + 0At° o i) 



c) 


FUx,Q*)/a 


0.0014 < x < 0.1 


O.K x < 0.87 


OPAL 
LO 
NLO 


0.180 ± 0.059 ± 0.026 ± 0.012 
0.065 

070 + °- 011 (0 044 + °- 004 + 026 + a007 ) 
u - u ' u -0.008 l u - u ^- 0.003 ^ U - UZD - 0.005/ 


0.084 ± 0.028 ± 0.010 ± 0.004 
0.084 

0.099 tofi (0.09818:8^ + 0.001) 



Table 1: a) The cross-section a D * in the restricted region, b) the cross-section <r(e + e~ — > 
e + e~ccX), extrapolated using Monte Carlo, and c) the structure function F^ c at Q 2 = 
20 GeV 2 . For the measurements the quoted errors are due, respectively, to statistics, sys- 
tematics and extrapolation. The statistical correlation between the two data points is —28%. 
The measurements are compared to several theoretical predictions. The HERWIG predic- 
tion in a) uses the value of /(c — > D*) contained in the program. The numbers in brackets 
are separate values for the point-like and hadron-like components, respectively. See text for 
further details. 
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Am [GeV] 



Figure 2: Distribution of the difference, Am, between the D* and D° candidate masses for 
both decay channels combined. The data are shown as points with statistical errors. The 
histogram represents the combinatorial background estimated using events with a wrong- 
charge pion for the D* — > D°7r decays. The curve is the result of the fit to a Gaussian signal 
plus a power-law background function, as explained in the text. 
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Figure 3: Background subtracted data distributions compared to predictions based on the 
HERWIG and PYTHIA programs. For HERWIG the hadron-like component is shown as 
a shaded histogram on top of the open histogram representing the point-like component. 
The Monte Carlo predictions are normalised to the data luminosity. Shown are a) Q 2 , b) 
the visible invariant mass, W vis , c) x vis , d) the charged multiplicity, A^, e) the transverse 
momentum of the D* with respect to the beam axis, pf*, and f) the absolute value of 
the pseudorapidity |^ D |. The data are shown by points with error bars representing the 
statistical error. 
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Figure 4: OPAL results for a) the cross-section cr(e + e~ — > e + e~ccX), with 5 < Q 2 < 
100 GeV 2 and b) for the charm structure function of the photon divided by the fine structure 
constant, F^x, Q 2 )/a, at Q 2 = 20 GeV 2 . The data points are the results obtained with 
the HERWIG Monte Carlo model. The outer error bar is the total error and the inner error 
bar the statistical error. The data points in b) are placed at those x values that correspond 
to the average predicted within a bin. The data are compared to the calculation of 
performed in LO and NLO. The band for the NLO calculation indicates the theoretical 
error from uncertainties in the charm quark mass and renormalisation and factorisation 
scales. In a) the cross-section prediction of the HERWIG Monte Carlo model is also given, 
b) also shows the prediction of the GRS-LO parametrisation for the structure function at 
Q 2 = 20 GeV 2 and its point-like component separately. 
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